Universal description for different types of polarization radiation 
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When a charged particle moves nearby a spatially inhomogeneous condensed medium or inside it, 
different types of radiation may arise: Diffraction radiation (DR), Smith-Purcell radiation (SPR), 
Transition radiation (TR), Cherenkov radiation (CR) etc. Along with transverse waves of radiation, 
the charged particle may also generate longitudinal oscillations. We show that all these phenomena 
may be described via quite simple and universal approach, where the source of the field is the 
polarization current density induced inside the medium by external field of the particle, that is 
direct proof of the physical equivalence of all these radiation processes. Exact solution for one of 
the basic radiation problems is found with this method: emission of a particle passing through a 
cylindrical channel in a screen of arbitrary width and permittivity e(u>) = e' + is". Depending on 
geometry, the formula for radiated energy obtained describes different types of polarization radiation: 
DR, TR and CR. The particular case of radiation produced by the particle crossing axially the sharp 
boundary between vacuum and a plasma cylinder of finite radius is also considered. The problem 
of SPR generated when the particle moves nearby a set of thin rectangular strips (grating) is solved 
for the arbitrary value of the grating's permittivity. An exact solution of Maxwell's equations for 
the fields of polarization current density suitable at the arbitrary distances (including the so-called 
pre-wave zone) is presented. This solution is shown to describe transverse fields of polarization 
radiation and the longitudinal fields connected with the zeros of permittivity. 

PACS numbers: 41.60.-m; 29.27.-a; 52.59.-f; 78.70.-g 



I. INTRODUCTION 

A charged particle moving nearby a spatially inhomo- 
geneous condensed medium (target) or inside it may pro- 
duce different types of radiation: Diffraction radiation, 
Smith-Purcell radiation, Cherenkov radiation, Transition 
radiation, Parametric X-ray radiation (PXR) etc. Along 
with transverse waves of radiation, the particle may also 
generate longitudinal waves. From the macroscopic point 
of view all these phenomena arise due to dynamic polar- 
ization of a medium by the external field of the moving 
charge and may be considered as the manifestation of the 
phenomenon known in the microscopic theory as polar- 
ization radiation (PR) [l[. Small radiation losses due to 
PR determine the wide range of its applications: from the 
new techniques of non-invasive beam diagnostics to 
the new sources of THz radiation 

However, up to now there is no universal theoretical 
description of all these radiation processes. The classi- 
cal Ginzburg-Frank solution for TR is known to include 
CR contribution [1 ll ] , but the corresponding solutions for 
DR and SPR generated on the surfaces of a finite permit- 
tivity (and therefore including possible CR) are still ab- 
sent, even for the simplest targets. Available results are 
bounded by TR and DR in X-ray range [H, [l3| , where 
the value of permittivity is close to unit: e(ui) — 1 <C 1, 
and the difference between various types of PR practi- 
cally vanishes (for example, between CR and PXR 14]). 
Moreover, in a real experiment one has different types of 
PR generated simultaneously. The absence of possibil- 
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ity to estimate, for example, a contribution of DR to CR 
generated on the targets of finite dimensions hampers the 
use of them in beam diagnostics [3, [l5| . The problem of 
coherent TR produced by the charged particles at the 
boundary of plasma occupying the finite space is of im- 
portance for development of new types of THz radiation 
sources [t| [l(| ■ But most available solutions are based on 
the presentation of the radiation from a sharp boundary 
between vacuum and plasma as from the "vacuum - ideal 
conductor" boundary (see, for example, flU [l7jV Such 
a model is rather approximate, and a more rigorous con- 
sideration is required. Finally, some recent publications 
were devoted to demonstration of an analogy between the 
theory of TR and DR and the one of the wake fields in 
accelerators [J, Il8l420l | . But there is still no good theory 
describing all the types of polarization radiation and the 
wake fields generated in a target of a given shape and 
finite permittivity s(uj) = s' + is". 

The exact solution of Maxwell's equations describing 
all the types of radiation mentioned as polarization ra- 
diation, is necessary not only for the numerous applica- 
tions of PR in accelerators, beam and plasma physics, 
but for more clear understanding of the general physi- 
cal properties of PR itself. For instance, TR is usually 
studied solving homogeneous Maxwell's eqs. with ap- 
propriate boundary conditions [ill H3- Energy losses 
on Cherenkov radiation in an unbounded medium is ob- 
tained as a part of the total polarization losses cor- 
responding to the excitation of transverse oscillations 
[2ll [22| . In most cases, CR and TR are considered 
separately from each other (though there exist several 
works directly devoted to the investigation of interfer- 
ence between them), but several authors suppose that 
the mechanism of their generation is the same (unfortu- 
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nately, without an adequate proof of this suggestion, see 
e.g. [23|). On the other hand, DR and SPR generated on 
the ideally conducting surfaces are often represented as 
the field of the surface current induced by external field 
of a particle (by analogy with diffraction of the ordinary 
plane waves, see e.g. j24l427l |). Though this formulation 
seems to be correct, it was recently shown that diffraction 
of the particle's field (not the one of a plane- wave) cannot 
be described with the use of a surface current with only 
two tangential components [28]. It turned out, that ex- 
act expression for the surface current density induced by 
the particle's field on an ideally conducting surface has 
all three components including the one perpendicular to 
the screen (generalized surface current |28j|). 

Due to the fact that such different ways are used to 
obtain the formulas for different types of the same po- 
larization radiation, the physical origin of these types 
of radiation is often just not indicated. In general, it 
is not enough to say that CR arises when the particle 
velocity is larger than that of light in this medium, as 
many authors still do. It is necessary to indicate the 
source of radiation, that is polarization current density 
in this case (see below). The well-known Cherenkov re- 
lation between angle of emission, particle velocity and 
permittivity is just a condition of the constructive in- 
terference of waves emitted at each point of the particle 
trajectory. The problem of the physical equivalence (or 
non-equivalence) of TR and CR eventually is connected 
with the fact, that exactly the source of radiation is often 
not indicated (because these are homogeneous equations 
that are solved in the case of TR). In what follows, we 
shall present the source of radiation leading directly to 
the well-known formulas for TR intensity including possi- 
ble CR contribution, that is direct proof of their physical 
equivalence. The term "equivalence" here should be un- 
derstood in the way, that both TR and CR, as well as 
DR, SPR and the longitudinal fields, physically arise due 
to polarization of a medium by the field of the moving 
particle. Therefore, the difference between these radi- 
ation processes is completely "kinematic" , and just the 
presence of interface of some shape leads to the changes 
in the radiation characteristics. 

The organization of the paper is as follows: in SecliTl 
we find the exact solution of Maxwell's equation for mag- 
netic field of polarization current density induced inside 
a target of a given shape and permittivity by the field of 
external source and the field of PR itself. In Sec lIIII we 
consider the problem of radiation generated by the par- 
ticle moving through cylindrical channel inside the layer 
of finite width and permittivity. It turns out that the 
energy radiated through the plane interface into vacuum 
may be easily found using the well-known reciprocity the- 
orem [21|. In Sec II VI we find the solution for the problem 
of polarization radiation produced by a particle cross- 
ing the boundary between vacuum and a plasma cylin- 
der of finite radius. Sec|V] is devoted to solution of the 
Smith-Purcell radiation problem arising while the par- 
ticle moves nearby a thin grating of finite permittivity. 



Finally, in Appendix the exact expressions for the fields 
of polarization current density applicable at the arbitrary 
distances (including the so-called pre-wave zone) are de- 
rived and analyzed. 



II. POLARIZATION RADIATION FIELD 
INSIDE A MEDIUM 

When a charged particle moves nearby a condensed 
medium or inside it, its own field polarizes atoms and 
molecules, and the secondary field known as polariza- 
tion radiation arises. In the microscopic theory this phe- 
nomenon is often called "atomic bremsstrahlung" (see 
e.g. [l[ and references therein). This secondary field 
propagating inside the medium in the form of transverse 
wave or exciton also polarizes atoms, that is why the 
source of PR is the polarization current density linearly 
depending upon the total field including external field of 
the movin g particle E° and the field W° l of PR itself (the 
medium is non- magnetic throughout the paper: [i = I; 
Gaussian system of units is used): 

j pol =a(E + EP° l (j pol )). (I) 

Note that this is the conventional expression for the 
macroscopic current density induced inside a medium 
used when averaging microscopic Maxwell's eqs. (see 
e -g- [HI H3)- Here conductivity is a given function of 
permittivity possessing frequency dispersion: 

o-(r,u) = ^-.(e(r,u)-l). (2) 

For non- magnetic media, Eq.([T]) just means that the 
field of polarization radiation is generated by electric 
dipole moments with effective volume density P(r,u;) = 
j pol (r,u)/(-iw) = (e - 1)(E° + W o1 )/{Att). 

Equation (fT|) is obviously integral, and neglect of the 
integral term E po/ (j po ;) is possible only for X-ray radi- 
ation, when the secondary re-emissions may be ignored 
due to smallness of parameter e(w) — I <^ I [H,[l3|- How- 
ever we can get rid of this term replacing it to the left 
side of "vacuum" Maxwell's equations, that leads to the 
replacement uj/c — > ^/e{t,ijj)lu/c. Maxwell's equations in 
this case may be reduced to the form: 

(A + e(r,w)^)HP°'(r,a;) = -— (a(r,u;)curlE 

-(E +E^) x grada(r,w)), (3) 

which still contains W° l in the right side. 

If the spatial dependence of conductivity <r(r,a;) is 
step-like: cr(r, uj) = 0(z)cr(w) (that corresponds to the 
plane interface between vacuum and medium) or has the 
more complicated form corresponding to a target of a 
given shape, it is easy to see that the last term in the 
right side of © may be expressed through the boundary 
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volume of the medium (target) Vt- An exact expression 
for electric field strength and its analysis is provided in 
Appendix. 

Note that this expression is the exact solution of 
Maxwell's eqs. with current density (JT]) in the right side. 
The integral term in the right-hand side of (TTJ resulted 
only in appearance of y/e(oj) in exponent. Physically, 
such procedure means "renormalization" of the field act- 
ing on each atom and molecule because of the presence 
of other atoms (that is macroscopically characterized by 
the permittivity e(lu)) and leads to the replacement 



Figure 1: Generation of PR by a charged particle moving 
through a channel in the layer 
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E°, uj/c -> VeMw/c, 
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conditions on the interface. For example, in the former 
case (n is a unit normal to the interface): 

(E° + E poi ) x grada(r, u) = a(u)5(z)E° x n, (4) 

i.e. tangential components of electric field are continuous 
on an interface. Therefore solution of Eq.Q may be 
written straightforward: 



H pol (r,u) =curl- / a(uj)E°(r', u) 



-dV,(5) 



where the region of integration is reduced only to the 



as formula ([5]) shows. The last is correct for any shape of 
the target, and in spite of its rather simple form describes 
all the types of PR inside a medium (see below) . 

To demonstrate how this approach works, we con- 
sider one of the basic radiation problems: emission of 
a charged particle moving uniformly along the axis of 
the round cylindrical channel with radius a in a screen 
of finite thickness d = d\ + 0I2 and arbitrary permittivity 
e(uj) = e' + is" (see FigfT]). In this geometry two types of 
PR arise: DR and CR, but if the hole radius approaches 
zero DR ought to transform into TR. Let us check these 
considerations. According to (J5)) PR field in the wave 
zone is found as (vector H R denotes the far-field part of 
H poi , see Appendix): 
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Here we use conventional expression for the field of uni- 
formly moving charge e with energy 7 = 1/ \J\ — (3 2 [29| 
(our coefficient in definition of Fourier transform differs 
from that in [29[ , but the final result for radiated energy 
does not depend on it of course). In (J7J also: K 0.1 are Mc- 
Donald functions, r' = {p',z'} — {p' cos (j)' , p' sin </)' , z'} 
is the vector of volume integration, v = {0, 0, v} is the 
particle velocity vector, and the wave vector in medium 



is k = £y/e(u)u}/c, where 

r = r{sin0cos</>,sin0sin</>, cos 8}. (8) 

Integration in ([7]) is performed using well-known pre- 
sentations for cylindrical functions [13] (see also such in- 
tegrals in Ref.[31|), and the final result for PR field is: 



R , , euJe(e -l)ae ir ^ u / c e-' ldl ^ -v^cose) _ e id 2 ^ (p -^cose) 

H (r, oj) = s ; — ; = — ^ ; — = -777- \ sin 6, — cos 6. \ 

V ' ' 2ttvcj 2 r (/3- 1 - ^cos9)(l - (3 2 + (/V^sm©) 2 ) 1 ^ Vl S 

x ( sin 6(7" 1 - /3-yy/e cos 6) J Q (a^-y/esm Qj K x (a^~) ~ (cos 6 + sin 2 6) Ji (a^ y/e sin e) K (a^~) ) , (9) 

I 



where J 0j i are the Bessel functions. 



First of all, we check the case of homogeneous medium: 
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d\,d,2 — > oo. Interference term is transformed as: 

oo 

e ^^f _ e id 2 ^f f e ifz dz = 2nif5(f). (10) 



As / = /3 1 — v / £~cos9, the zero of delta- function gives 
the ordinary Cherenkov condition: 



cos 9 = 



1 



(11) 



and the formula (|9|) describes CR generated in a channel 
in the unbounded medium 32|. 

In the limit a — > the following relations are fulfilled 
(arguments of Bessel functions are omitted) [30| : 



«7 

ciJqKi — > — , aJiKo — > 0, 
ui 



(12) 



and the energy radiated into unit frequency interval is 
found using @ as (the medium is transparent): 



^=/f£|H«|»«tt=*^- 

duj J y/e 7TC 



j sin 2 I3y/E cose) 

(l-/Vecos9) 2 



(e-l)(l-^ 2 -/Vecos9) 



l-/3 2 + (/V^sin9) 5 



sin 3 9rf9. (13) 



If the Cherenkov condition is fulfilled, the term in the 
square brackets is equal to unit, and the formula (|13p 
is reduced to the one known in the Tamm theory (see 
e.g. [33|]). For the layer of large width, interference term 
transforms into delta-function: 

sin 2 U%(l-p</ecaaQi) 

(l-/Vecos9) 2 T 

_> i-n 5(l- p^cosO). (14) 
2 v \ ) 

After this transformation integration in (1131) gives 



1 dW e 2 
d dio 



r w 1 



P 2 e(w) 



(15) 



that is the ordinary formula for CR spectral density in 
the unbounded medium (see e.g. [HI). 



III. THE ENERGY EMITTED INTO VACUUM 

Formula (J9j) has been derived for polarization radiation 
field taken inside a medium. In order to find radiation 
intensity outside the target (in vacuum), we cannot use 
the ordinary Fresnel laws of refraction because for ab- 
sorbing medium the radiating dipoles are concentrated 
in the vicinity of interface and the radiation field in this 
region does not correspond to the far-field (it is always 
so if e" ~ e'). To solve this problem, it is useful to ap- 
ply the so-called reciprocity theorem [2l[. The idea is 



in solving the inverse problem, we are to find the field 
refracted into medium if the field of the wave incident 
from vacuum onto xOy plane is determined by Eq.(|9]). 
In this case it is possible to use the ordinary Fresnel laws 
of refraction, and finding the solution for such inverse 
problem we can solve the initial problem, i.e. obtain the 
field refracted from medium into vacuum. Thanks to az- 
imuthal symmetry, we find the simple relation (in more 
detail see Ref.[H): 
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where 



sin (9 = y'isinS, 



(16) 



(17) 



9 is the polar angle in vacuum [21(, and fn is the Fres- 
nel coefficient. Note that no additional approximations 
were made when deriving Eq. (ll6[) . It should be also em- 
phasized that though the use of a polar angle in medium 
9 has a sense only for transparent media, all the results 
below have no restrictions on the imaginary part of per- 
mittivity. 

Thus in order to find the energy radiated in a unit 
solid angle and frequency interval according to relation 
d 2 W/dujdQ = cr 2 |E^ QC | 2 , it is enough to express formula 
© through the vacuum angle 9. Doing this and suppos- 
ing d\ = d, di = we get finally: 
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The expression derived describes all the types of PR gen- 
erated in a layer and emitted into vacuum: CR, DR or 
TR (if a = Q). 

In the limiting case of semi-infinite continuous medium 



(a — > and d — ► oo) the formula (fT5)) is reduced to the 
well-known Ginzburg-Frank expression for forward TR 
ELE1: 
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Here (as well as in (|18p) the pole /3\/ £ — sin 2 6* = 1 cor- 
responds to CR cone refracted into vacuum, but for the 
target of finite width this pole is removable (the medium 
is transparent): 
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and the radiation intensity is proportional to the width 
squared, as in Tamm theory of CR without channel (see 
e-g- [13, H3| and references therein). Note that integra- 
tion over angles leads to the linear dependence on d. 

The argument of sin 2 (d/d/) in (|2"01 is related to DR 
(TR) formation length in the weakly absorbing medium: 



A 



ttI/3- 1 - Ve 



(21) 



becoming infinity for the "pure" CR (when e" — 0). If 
the target thickness d is much less than df, we get for 
TR the same oc d 2 dependence even if Cherenkov con- 
dition is not fulfilled. If the opposite condition holds 
true d ^> df (thick target), the interference term again 
transforms into delta- function according to (|14l) . deter- 
mining the Cherenkov condition via the vacuum angle 
9 (because in this case CR is observed in vacuum, but 
not in the medium as in the previous paragraph). Far 
from the Cherenkov angle, sin 2 (d/d/) for the thick tar- 
get is highly oscillating and may be replaced by its mean 



value 1/2. So, in contrast to CR the intensity of tran- 
sition radiation (diffraction radiation) in a transparent 
medium almost does not depend on the target thickness 
when condition d^> df is fulfilled, see Fig. 2. 

In another limiting case of conducting target with 
e" » e' , exponential term in (Ti"%)) , ([2"U|) becomes damped, 
and dependence upon the target width disappears if con- 
dition d 3> A/-v/jef is fulfilled {skin- effect). In the case of 




20 

0(deg) 



Figure 2: Angular distributions of the total polarization ra- 
diation in vacuum. Parameters: 7 = 10, A = 1mm, a = 1mm, 
e' = 1.3, e" = 1(T 3 . 



ideal conductivity (e" 
forward DR: 



formula (fig)) describes only 
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This is exact expression for DR generated by a par- ticle moving through a round aperture in an ideally- 
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Figure 3: Angular distributions of the total polarization radiation in vacuum. Parameters: j = 10, A = 1mm, d = 40mm; a.) 
a = 5mm, e' = 1.3; b.) e = 1.3-M0.05 (dashed curve is multiplied by the factor 60). The angle of Cherenkov radiation reflected 



into vacuum is determined from the condition /3y e — sin 2 9 = 1. 



conducting screen, that was found with another methods 
e.g. in Refs. [U SI. 

Note that formula (fT5|) describes only forward po- 
larization radiation. The corresponding expression for 
backward emission is found by the substitution cos — > 
— cos 9 in © and applying the same refraction rule (TTBT) , 
see also [34|. 

In the general case, when the particle moves through 
the channel in the screen with arbitrary permittivity 
e(ui) — e' + is", formula (fT5|) completely describes this 
process, i.e. represents the total radiation intensity in- 
cluding Cherenkov and diffraction mechanisms. It should 
be noted that we consider refraction of polarization radia- 
tion from medium into vacuum through the flat interface 
laying in xOy plane (see Fig.l), so we do not take into ac- 
count the secondary reflections of the emitted radiation 
inside the channel, that seems to be relevant when the 
screen width d is large enough. However this approxima- 
tion is in fact insignificant, because Cherenkov radiation 
as is well known is not radiated through the interface 
parallel to the particle trajectory (see below). In other 
words, the total internal reflection takes place, and in this 
case CR is emitted into vacuum only through the surface 
perpendicular to the particle trajectory (if CR condition 
sin 2 6 = \e - l//3 2 \ < 1 is fulfilled), see e.g. H3- 

The angular and other distributions of the energy 
losses due to different types of polarization radiation are 
shown in Figs l3l4l5l As may be seen in Figj3^i.), increase 
of imaginary part of permittivity leads to CR suppres- 
sion. Increase of the channel radius results in transform 
of TR into DR and in decrease of CR contribution, see 
Fig|3j3.). Fi^4] shows, that the presence of sin# in the 
argument of the oscillating Bessel functions in Eq. fTE)) 
leads to the different dependence on a for intensity of 
diffraction radiation (small 0) and the one of Cherenkov 
radiation (large 9). The difference between them reveals 
itself also in the dependence of the radiated energy on 
the Lorentz factor of the particle, see FigJSJ In particu- 
lar, this figure shows that for the large enough values of 
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Figure 4: Dependence of polarization radiation intensity on 
the hole radius. Parameters: 7 = 10, A = 1mm, d = 50mm, 
e' = 1.3, e" = 10~ 3 . 
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Figure 5: Dependence of polarization radiation intensity on 
the particle energy. Parameters: A = 1mm, a = 2mm, d = 
50mm, e' = 1.3, e" = 10" 3 . 



7 (actually, when 7A/(27r) 3> a), energy dependence of 
DR becomes quadratic as for TR [ill ]. 

It is clear from the figures that all three types of radi- 
ation (TR, DR and CR) have the same physical origin, 
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and the boundary between them is rather relative. Note 
that the physical equivalence of TR and DR was shown 
for an ideally conductin g ta rget in the ultrarelativistic 
limit relatively recently [35j |. Here we proved this for 
the case of arbitrary permittivity and arbitrary particle 
energy that allowed to consider possible contribution of 
Cherenkov radiation. 



electrons (125) . zero damping and heavy ions: rrii 3> m e , 
the last relation is reduced to ([M]) under the condition 
of small density 



or high frequency 



n < 



(27) 



IV. POLARIZATION RADIATION FROM A 
PLASMA CYLINDER 

Study of radiation generated by the charged particles 
crossing the boundary between a channel of the dense 
plasma and vacuum is of importance for development 
of new THz radiation sources [To| | . TR from such a 

17| as the one pro- 



1G 



boundary was modeled in Rcfs. 
duced at the boundary between vacuum and an ideal 
conductor. This was justified by the fact that the value 
of the electron plasma permittivity 



1 



u p(e) 



(23) 



tends to — oo when w p ( e ) 3> ui (dense plasma, see 
Ref.[ll). 

Actually, for frequency range below the optical band, 
this formula is not applicable for collisionless plasma, be- 
cause the effects of the spatial dispersion of permittivity 
are ignored. But even neglecting by the spatial disper- 
sion, this formula works only for small densities (w p ( e ) is 
of order or less than w), because it is based on the sup- 
position of the linear dependence of permittivity on the 
density of electrons n e : 



e(w) = 1 + 47rn e a e (o;), 



(24) 



where a e (uj) is the polarizability of electrons (electrons in 
plasma or conduction electrons in a metal, see e.g. prob- 
lem 11.17 in [11). In the simplest model of oscillators 
with zero own frequency (cold plasma) and attenuation 
factor r e , the polarizability of such electrons is (see e.g. 

[HH) 



a e ((jj) 



1 



(25) 



ifr e 

notation uj 2 , N 

p(e) 



m e u 2 + iT e uj ' 

(i.e. neglecting by Landau damping), and with 
47m e e 2 /m e , the formula (|2~4"|) coincides 

with flU). ' 

The more rigorous study allowing to consider high den- 
sities is based on Clausius-Mossotti formula for an elec- 
trically neutral isotropic medium of weakly interacting 
particles: 



s(uj) 

Here na(u>) = 



1 + ^Lna(w) 



= 1 



(26) 



1 - fna(u) 1 - fna(w) 

^2n s a s (u>) represents the sum over all 



subsystems of the medium being studied. For electri- 
cally neutral (as a whole) plasma with polarizability of 



U >• U) p ( e ) , 



(28) 



where r e = e 2 /(m e c 2 ) is the classical electron radius. 
As the opposite case w <C WpU) does not led to (|24|) . it 
is impossible to apply the formula (1231) for high density 
plasma. 

Consider the partially ionized electrically neutral 
medium consisting of free electrons, bound electrons and 
ions. Polarizability of free electrons is determined by 
Eq. (l25"]) . The analogous expression for bound electrons 
takes into account the possible own frequency of the os- 
cillators uiq (only one for simplicity): 



a&. e .(w) 



1 



m e oj 2 + iTiUj 



(29) 



Since on oc m i , we shall neglect by the contribution of 
ions due to their big mass: m, 3> m e . So the formula 
(j26| yields: 



2 

1 - - 



J p(e) 



Mi) 



3 oj 2 + iT e ui 3 ui 2 + iTiU 



1 

1 + - 



J p(e) 



1 



,(30) 



Mi) 



3 oj 2 + iT e bj 3 w 2 + iTiUj 



where u 2 ^ = Airriie 2 / ra e is denoted, and rij is the density 
of ions with bound electrons. If the term uj 2 + iTiUj — a; 2 
does not turn into zero, and the density of the partially 
recombined ions m is much less than the one of free elec- 
trons n e , the expression for permittivity describes only 
electron plasma: 



e{u>) 



'Me) 



+ o W p(e) + 



(31) 



Here, the term kui 2 ,^ in the denominator appears be- 
cause of the taking into consideration the so-called lo- 
cal field effect in Clausius-Mossotti formula. In fact, for 
high density of a medium (i.e. when condition (|27p is not 
fulfilled) the field acting on each atom or an electron in 
plasma does not coincide with the commonly used macro- 
scopic field E = E° + W° l , and dependence of permittiv- 
ity on the density becomes non- linear as Eq. p5|) shows 
(see also [22],|36j]). In this case, even the formal transition 
w p (e) — > oo in (|3ip (having no physical sense whatsoever, 
because the formulas (|2~4"|) . (121)1) were derived for a set of 
weakly interacting particles) results in the limit value of 
permittivity s(uj) — > —2, but not — oo. It means that the 
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modeling of high density plasma as the ideal conductor 
used in Refs.[l6l. Il7| appears as the rather rude approx- 
imation. On the other hand, for the special problem of 
PR generated by a relativistic particle in a plasma cylin- 
der, the difference between the radiation characteristics 
for the value of permittivity e ~ —1.5 and for the case 
of ideal conductor (e = ioo) is significant only for large 
angles of emission 8 ^> 7 -1 (see below). 



Consider an electron crossing axially the sharp bound- 
ary between the plasma cylinder of radius a and vacuum. 
The forward transition radiation may undergo diffraction 
due to the finite radius of cylinder. Such a problem is 
"complementary" to the one has been solved in the pre- 
vious Sections. So, the radiation field inside the cylinder 
is found using Eq.© (integrating over p' from to a) as: 



fZ( p _ -\\ ir^/su) /c -id-i_^{p - ^/e cos 0) _ id 2 ~ (ft -^/ecosQ) 

H fi (r,u,) = H* fl - © = 6V T lj - = * , ta r . ^ a , {sM,-co B &0} 



2ttc r (/3- 1 -Vecos8)(l-^ 2 -(-(/3V£sine) 2 ) 

f sin 6(1 - P 2 - ^VecosO) - a— sin 6(1 - /3 2 - /V^ cos 0) J fa- VesinOW fa—) + 
V ^7 V c / V vy/ 

+a-^ sin 6(cos 6 + fi^fe sin 2 6) Ji (a^ ^/i sin 9) if ( a ^~) ) » ( 32 ) 



Then, following the procedure completely analogous to 
the one used in Sect lllli we consider refraction of the 
emitted radiation through the plane cylinder "top" lay- 
ing in xOy plane in Fig[TJ that restricts us in what fol- 



lows within the radiation angles 9 <C w/2. Applying the 
reciprocity theorem, we come to the following formula for 
the energy radiated into vacuum (semi-space z > 0) from 
cylinder of length d (d\ = d, d% = 0): 



2 P 2 cos 2 6 -tdf f^-Ve-sm 2 e) 
- — — — — — e V / — 1 



d 2 W 



ckudfl 7r 2 c (1 — P 2 cos 



£ - 1 



(1 — flye — sin 2 9){e cos 9 + \/ e — sin 2 



f sin 0(1 - /3 2 - /3a/e- sin 2 0) - a— sin (9(1 - /3 2 - pVe- sin 2 6>)J fa- sin 6*") if x fa— 
V ^7 V c / \ 1)7 

+a-^(V £ - sin 2 9 + /3sin 2 9)J 1 (a— sin flWofa— 

vp V c /V W7 



(33) 



It is not difficult to see, that the limit a — > 00, d — > 00 
leads again to the ordinary Ginzburg-Frank formula for 
forward TR (with contribution of CR) at the infinite 
boundary (|19p . Substitution of formulas for plasma per- 
mittivity $Z3§ or (|31[) into (|33l) makes the last expression 
describing PR generated at the boundary of vacuum and 
cylinder of cold electron plasma. For the high density 
plasma and THz region of frequencies, one should use 
formula (|31l) . For the high frequencies 10 ^> uJ p ( e ) (ultra- 
violet and X-ray regions), it is possible to use the simplest 
expression (|23|) . see also [l3j]. We emphasize, however, 
that for THz range formula (|31[) may be considered as 
just the first correction to (|23p taking into account the 
possible local field effect in the high density plasma, but 
not taking into account the spatial dispersion. 

The angular distributions of the radiated energy are 
shown in FigJBJ The peak of TR is shifted to the large an- 
gles region due to small radius of cylinder a <C 7A/(2tt). 
The negative value of permittivity results in the absence 



of CR, as the term 1 — /3y/e — sin 2 9 does not turn into 
zero. This feature is known to be the one of the radiation 
of an electron moving in the cold isotropic plasma, where 
only the longitudinal waves may be emitted 37]. As is 
clear from FigJS]b.), the significant difference between the 
radiation characteristics for plasma with e = —1.5 and 
for the ideal conductor has a place in the region of large 
angles of emission 9 3> 7 _1 . But for the angles 9 ~ tt/2 
the model, strictly speaking, is not applicable. 

If the radiation properties at large angles 9 ~ 7r/2 are 
of interest (for instance, in the non-relativistic case), one 
should consider refraction of the exact field (|32|) through 
the concave surface of cylinder (not through the back 
surface, like in (1331) ). The last problem may be simplified 
in the two limit cases: 1.) a ^ A and 2.) a ^> A. In 
the first case, the channel represents the small spatial 
inhomogeneity (see also the next Section), and applying 
of reciprocity theorem is not necessary. In the second 
case, the surface of cylinder is almost fiat for such waves, 
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and refraction through this surface may be calculated as 
through the flat one. 

It is worth noting the special case of partially ionized 
gas with some quantity of bound electrons possessing 
their own frequency ujq. If the damping constant I\ in 
([50)1 is negligibly small, and the frequency of interest w 
is close to wo, the terms of bound electrons in (|50j) dom- 
inate: 

e(u) « 1 ^ -> -2, (34) 

when to — > wo (though neglect by the local field effect 
would led again to -co). Such a model allows to con- 
sider PR at the frequencies which are close to the natural 
ones of the medium. Note that DR on such a resonance 
frequency was considered for the first time in Ref . 38] . 

V. SMITH-PURCELL RADIATION AS 
POLARIZATION RADIATION 

Another type of polarization radiation is generated 
when a particle moves uniformly at the distance h par- 
allel to a periodic surface (grating) 39]. The simplest 
target for this type of radiation is a set of N rectangular 
oo x b x (d — a) strips made of material with permittiv- 
ity e(uj) and separated by vacuum gaps. Here d is the 
grating period, a is the vacuum gap width, see Fig[7] 
Physically, there is no difference between Smith-Purcell 
radiation and the ordinary diffraction radiation. Fur- 
thermore, they both (as well as TR, CR and PXR) arise 
because of polarization of the medium under the action 
of external field of the particle. And only periodicity of 
the grating changes the radiation spectrum, that makes 
SPR different from the ordinary DR. 

Equality ([5]) represents the exact solution of Maxwell's 
equations describing all the types of PR inside a medium. 
But it was deduced only for a "simply connected do- 
main" , and in the general case it is not correct for the 



completely separated regions. It would be correct for 
a grating with continuously deformed profile, for exam- 
ple, rectangular or sinusoidal one. However, if the gap's 
width a is much less than a wavelength in the transpar- 
ent medium a <C A/ yj £(oj), the gaps reveal themselves 
as the small spatial inhomogeneities (see also Ref. (4p|j ) . 
The scattering cross-section of a plane wave on such "par- 
ticles" is proportional to their polarizability [2l|, which 
is obviously equal to zero as e vac = 1. As regarding a 
conducting grating, in this case the waves are emitted 
only from the thin layer near the surface of a strip (skin 
depth) , and they do not reach the interface as the trans- 
verse plane waves. Therefore, condition of applicability 
for Eq.© becomes just a -C A. 

In order to find the energy radiated into vacuum, one 
should consider refraction of a spherical wave on the 
rectangular surface of a strip. The last problem may 
be simplified supposing that condition b -C (d — a) is 
fulfilled, that excludes from the following consideration 
the angles of radiation which are too close to the sur- 
face (6 — > 0, 7T, see similar calculations for a rectangular 
screen in Ref. [34] ) . In this case, it is possible to consider 
only refraction on the upper (or lower) plane of a strip. 
Applying again the reciprocity theorem, one should find 
the field refracted from vacuum into medium through the 
surface of a strip. But according to the Smith-Purcell re- 
lation (see below) , the last one has a dimension of order of 
a wavelength (d ~ A), and therefore cannot be considered 
as the homogeneous flat surface. So, the only case allow- 
ing to apply the ordinary Fresnel laws of refraction is the 
geometry when the same condition a <C A is fulfilled, and 
the surface of the grating is almost homogeneous for such 
waves (see e.g. [4t|). However, this condition does not 
bound the values of permittivity e{uS) = e' + is" being 
considered, so it is possible to develop the exact theory 
of Smith-Purcell radiation in the sense of taking into ac- 
count the real properties of the grating material. 

The radiation field in the wave zone is found from 
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y 



Figure 7: Generation of polarization radiation by a charged 
particle moving nearby a grating. 



Eq.([5]) in the following way 



Li.' 



H> >W ) = -( £ ( W )-1). 

b 



ir-\J e{uj)uj / c 



kx / dz' J dy' E°{k x ,y',z',Lj)e- 



,(35) 



where the integral over z' represents a sum of the inte- 
grals over each strip of the grating. The Fourier trans- 
form of the particle field now has the form (see e.g. (3lj |) 



E u (k x ,y',z',Lj) 



-le e 



2nv ^i + e^ej- 



^ v /l + e(/37e x )^ 1 ~ 1 } e - (y ' +h) ^^ 1+£( ^ e * )2 ,(36) 

Integration over z' gives the well-known product of inter- 
ference terms: 



dz 



2sin(^(^ - k z ) 
- - k z 

V * 

sin 



3in(i(^-fe,)) 



FstripFN- (37) 



Substitution of (|36| and f|37j) into (|35|) leads to the fol- 
lowing expression for radiation field (dependence of e on 
lo below is implied): 



2 ir^/e — 

exp ( - b(i^e y + ^^l + e{he x ) 2 )) - 1 sin (^(H - fc,)) sin (iVf (| - fc z ) 



e-^V^+e^-r^) 7 , (38) 



If the particle moves above a homogeneous medium 
(N = 1, <i— a — )• oo) the interference term is transformed 
into delta-function: 



sin -*,)) 



(39) 



determining the Cherenkov condition. However, this field 
propagates inside the medium and is not radiated into 
vacuum through the upper (lower) surface (see below and 
also jal). 

The next step is to find the energy radiated into vac- 
uum using the reciprocity theorem. As we already used 
the condition a -C A, one may consider the grating as 
the homogeneous interface and apply the same rules of 
refraction as in the previous paragraphs. As the prob- 
lem being solved has no azimuthal symmetry, the final 
expression for the absolute value of the radiation field in 
vacuum becomes (see also 34] ) : 



IE, 



| 2 = Ti|#ff 



■T\\(\H, 



R\2 

y I 



<| 2 ), (40) 



where 



H ± = H*cost 



H? sin ( 



ff,f = tf^sinc 



iff cos. 



(41) 



are the components of magnetic field (|38[) perpendicular 
and parallel to the plane of incidence of the wave on the 
interface of the grating, and also: 



'Hi = 



Ih 



e 
Je 



2 sin 8 cos < 



e sin f cos < 



• + y/e- l + (sin< 



cos < 



2 sin 8 cos < 



y/e(sm 8 cos <j> + \J e — 1 + (sin 8 cos < 



(42) 



are the transmission coefficients expressed through the 
vacuum angular variables used in the Smith-Purcell ra- 
diation problems: 



eg = {sin 8 sin 0, — sin 8 cos 4>, cos 8} , 



(43) 



and fn,fE are the Fresnel coefficients [2l[. As may be 
seen in Fig|6l the azimuth angle <f> is counted from the 
negative direction of y-axis, because we consider radia- 
tion emitted into the upper semi-space y < 0. The az- 
imuth angle 4> in formula (|4"T]) belongs to the plane of 
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interface: xOz. Connection of this angle with variables 
(|4*3"| is found by equating the components of (|4"3")l and 
the ones written through the angles used in a problem of 
refraction (see e.g. [21| ) . The final result is: 



cos 9 



sin < 



i,COS< 



v/l — (sin 9 cos <f) 2 



The last step is to express the unit vector of radiation 
in the medium e = r/r in Eq. (|38[) through the vacuum 
variables according to the Snell's law 21] : 



{sin 9 sin <ft, — yjs — 1 + (sin 9 cos </>) 2 , cos 0}. (45) 



The minus before the square root is chosen so that the 
wave (|38|) is damped deep into absorbing medium. Re- 
turning to the delta-function determining the condition 
of Cherenkov radiation ([55)1 . one may note that its ar- 
gument written in the vacuum variables is proportional 
to 1//3 — cos 9 and never turns into zero. It means the 
absence of CR transmitted through the surface paral- 
lel to the particle trajectory, as we already indicated in 

sectnm 



Now we can write down the final expression for the 
energy radiated in the upper semi-space: 



d 2 W e 2 


(/3sin#cos^) 2 sin2 (V7(!" C0S(? ) 


)sin 2 (iV|H(I 


-cos 0)) 


dudfl 7T 2 C (1 - f3 2 - 


V (/3sin<9sin0) 2 )(l - (sinflcos^) 2 ) (1 - /3cos6>) 2 


sin 2 


cos efj 



(e - l)(exp{-&^( v /l + (f3j sin 9 sin <j)) 2 - i^y/e - 1 + (sin 9 cos (j)) 2 )} - 1) 



a/1 + (/37sin6»sin0) 2 - ifliy/s - 1 + (sii 



(7 1 cos8 + /^(sin^sin^) 2 )^ — 1 + (sin 9 cos 4>) 2 + i(l - (sin cos 0) 2 )-\/l + (/37sin#sin< 



^2 2 



e sin f cos ( 



+ y/s-l+ (sii 



+(sin0sin^) 2 (7 x — ,87 cos ( 



, 1 - (sin cos </>) 2 + I ^/e — 1 + (sir 



^212 



\y/e(sin9 coscj) + \J e — 1 + (sin6>cos0) 2 )| 2 



sin y sin < 



(46) 



If the number of strips is large enough (N 1) the 
interference term squared gives another delta-function: 



5in 2 (j\ffs(i _cos0) 



sm 



(^(1_ CO S0)) 

27riV ^ <5 (d- ( - - cos (?) - 2irm) . (47) 

™ — i c p 



Its zeros give the well-known Smith-Purcell relation [3£ 



d /l 

A m = — - - COS I 

m Vp 



(48) 



Although the formula (|46| is rather complicated, it 
permits to do the following factorization 



d 2 W d 2 W 



diudVl diudVl 



r strip r N ') 



(49) 



where d 2 W /r /dajdil| sc is the energy radiated while a charge 
moves above a semi-infinite screen. In other words, the 
first factor in (|4"!)f describes DR generated on a semi- 
infinite screen of the finite permittivity (conductivity) 
and height. The second factor describes interference of 



the waves emitted from the left and right sides of a strip. 
Finally, the third factor takes into account the interfer- 
ence of waves emitted from all N periods of the grating. 
This factorization is known in literature for the case of 
an ideally conducting grating [4l| . In the present theory, 
it appears as a result of condition a <C A used when de- 
riving Eq. (|4^|) . In fact, this limitation is quite general, 
because if the opposite condition o> A holds true, it is 
impossible to consider radiation from a single strip (even 
the ideally conducting one) to be entirely independent 
of radiation from the other one. In this case, the sec- 
ondary re-reflections of radiation may be of importance, 
and the formula for the radiated energy would not have 
the structure (|4"9")l . Note that such factorization does not 
have a place also in the van den Berg's model of SPR 
(calculations within this model for the grating being con- 
sidered were performed in Ref . (42j ) . It is so, because in 
this model the limitation b <C (d — a) was not used. Un- 
fortunately, it was done only for the ideally conducting 
gratings. 

To begin analysis of the formula (|4l)|). let us consider 
the simplest and the most investigated case of ideal con- 
ductivity: e" — > 00. Performing the corresponding tran- 
sition, we have 
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d 2 W e 2 ^- 2 (1~ (sin0 sine/)) 2 ) + 2/3 cosd(sm9sin(j)) 2 + (f3^) 2 sm 4 e sm 2 < 



1 + (,$7 sin (9 sin 0) 2 
sin 2 -cose)) sin 2 -cos 0)1 , ; 

(l-/3cos0) 2 sin 2 (f^(± -cos0) N 



where dependence upon the grating height b has disap- 
peared due to the skin-effect on the upper facet of each 
strip. For the large number of strips, we can integrate 
this expression over frequency using (|47p . The result for 
radiation on the m-th harmonics is 



1 dW m 

N dn 



2e 2 (3 
ird 



sin' 



i?( 7) < 



{I -/3 cose) 3 



- w l + ([3y sin 9 sir 



(51) 



where we use denotation R(jj, <fi) for angular part (the 
first line) of formula (|5*0| . 

It is essentially that the formula (|50p may be derived 
using completely different way. As we already noted, 
there exist several models of SPR for the ideally con- 
ducting gratings based on the presentation of SPR as a 
field of a surface current density induced on the grating 
by the field of the particle. But the common drawback of 
these models consists in the neglect of the normal to the 
interface component of the surface current (see details in 
28]). Without taking into account this component, it is 
impossible to obtain even the simplest Ginzburg-Frank 
formula for intensity of transition radiation on an ideally 
conducting plane. On the contrary, it is the exact ex- 
pression for the surface current density found in Rcf.]28] 



isurf 



C 
Z7T 



n x E° 



(52) 



that leads to the well-known results for TR (including 
geometry of oblique incidence). Here, n is the normal 



vector to the surface. This expression allows to obtain 
exact formulas for DR, and in the special case of the 
ideally conducting grating it leads to the very same for- 
mula (|50| for Smith-Purcell radiation (43|. Lastly, the 
part d 2 W/dujdQ\ sc of (jSTJ)) corresponding to DR from a 
semi-plane completely coincides with the one derived in 

Refs.gllllll- 

Comparison of SPR angular characteristics according 
to formula (|51l) and the ones within the known surface 
current model [26] without the normal component is pre- 
sented in FigJSl Note that analog of formula (f5"Tj) in the 
model 26} is obtained by the substitution: #(7, 6>, (f>) —> 
(/3sme) 2 . It is easy to see, that in addition to disappear- 
ance of the important angular terms corresponding to DR 
from a semi-plane, the extra /3 2 arises. Completely the 
same excess term also has a place in the formula for TR 
on an ideally conducting plane obtained with the method 
of the surface current without normal component similar 
to that of Ref.[26| (see in more detail [HI, Hl|). How- 
ever, in the ultrarelativistic limit (7 — > 00), the function 
,R(7, 9, <fi) in (|5ip tends to sin 2 e, and the minimum in the 
azimuthal dependence disappears (compare a.) and b.) 
in Fig J8]). Thus, for ultrarelativistic energies the use of 
the surface current model [26j is justified (for example, 
for 28.5 GeV electrons of the experiment [5j at SLAC). 
For moderate relativistic electrons (hundreds of MeV and 
less), one should use the exact expression (|5Tj) . Coinci- 
dence of both models in the ultrarelativistic limit is no 
surprise, because the model of a surface current with- 
out the normal component originates in the theory of 
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diffraction of the ordinary plane waves [H, [3l| . In the 
ultrarelativistic case, the field of an electron is almost 
transverse, that is why the methods of optics are quite 
applicable. We emphasize, however, the limitation of all 
the surface current models for the gratings with vacuum 



gaps: the gap width must be less than the wavelength 
studied. 

For the opposite case of transparent material of the 
grating (s" = 0), dependence on the grating height be- 
comes 



exp 



I try V 



'1 + (/^sintfsin^) 2 
b 



ifay/e - 1 + (sin 9 cos 4>)A j 



- 1 



= a( sinh 2 (- — Jl + (/3 7 sin6»sin0) 2 
V V 2 vi 



sin 



-J e — 1 + (sin 6 cos ( 

2 c 



(53) 



r 



For the values of b <C where sinh 2 (x) ~ x 2 <C 

1, this implies the almost periodical dependence of the 
radiation intensity as on the grating height b, as on the 
value of permittivity e = e' . However, for conducting 
material of the grating with e" 1, the radiated energy 
weakly depends upon the value of permittivity. 

In the general case with e(u) = e' + is", intensity of 
Smith-Purcell radiation for a conducting grating several 
times exceeds the one for the grating of transparent mate- 
rial, see FigjSJ On the contrary, intensity of the ordinary 
DR generated within small polar angles weakly depends 
on the value of permittivity. This forward diffraction ra- 
diation is generated on the grating as on a semi-infinite 
screen. 

It should be noted, that the well-known absence of ra- 
diation on even harmonics for a — d/2 predicted by the 
surface current models of SPR (see (fipTj) and (2(| |4l|) ap- 
pears just due to the fact, that all they are applicable 
only if condition a <C A ~ d is fulfilled. Argumentation 
presented in the beginning of this Section clearly demon- 
strates the origin of this restriction. In conclusion of 
the Section, we emphasize that solution for SPR derived 
above represents the exact solution of the problem un- 
der the conditions used: small vacuum gaps, thin strips. 



In order to find an analogous solution for a grating with 
continuously deformed profile (rectangular or sinusoidal 
one) it is not necessary to use the restriction of the small 
gaps. But in order to overcome condition b <C [d — a), 
we have to know the laws of refraction of a plane wave 
on a wedge having the right opening angle and made of 
material with e — e' + is" . Though we do not know the 
solution of such a problem, it may be useful to apply the 
reciprocity theorem at least for a conducting wedge using 
the well-known Sommerfeld's solution (see, for example, 

mi). 



VI. 



CONCLUSION 



As is shown in this paper, the well-known formulas for 
Cherenkov radiation, Transition radiation and Diffrac- 
tion radiation generated by the charged particles in the 
targets of the real dielectric properties may be obtained 
by solving Maxwell's equations with polarization current 
density |T]) in the right-hand side. Such universal de- 
scription proves the physical equivalence of the various 
types of polarization radiation. Moreover, this approach 
allows to calculate with a good accuracy characteristics 
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of TR, DR, SPR and CR generated in the targets of the 
given shape and finite permittivity e = e' + s", that is 
of relevance for development of new techniques of beam 
diagnostics and new types of radiation sources. In par- 
ticular, the solution found for SPR problem is applicable 
for the arbitrary energies of particles (in contrast, for ex- 
ample, to the model [26|) and the arbitrary value of the 
grating's permittivity. The model developed also allows 
to simulate the radiation characteristics from an electron 
beam transported through the special dielectric struc- 
ture, for instance, a cylindrical waveguide with a grating 
produced on the internal surface of the one. In this case, 
radiation will be emitted through external surface of the 
cylinder and may be detected without problems. Such 
a scheme provides minimal beam perturbation because 
of the azimuthal symmetry and low SPR energy losses, 
and it may be used for non-destructive beam diagnos- 
tics. Finally, the exact taking into account of the spatial 
dispersion of the target's permittivity may be also per- 
formed, that allows to develop more rigorous models, for 
example, for PR in the spatially inhomogeneous plasma. 

It should be mentioned that in order to detect CR in 
vacuum from a target with e' 3> 1 (i.e. when condition 
sin 2 9 = \s — 1 / /3 2 1 < 1 is not fulfilled), the normal vec- 
tor to the "exit" surface should coincide with CR direc- 
tion in the medium = arccos (l/(/3y / £))- Calculations 
for a target made of material possessing weak dispersion 
and having such inclined surface may be performed us- 
ing directly Eq.© and neglecting by refraction at the 
boundary. It allows to develop the new technique for 
the non-intercepting bunch length diagnostics based on 
the Cherenkov radiation generated in the targets of the 
special shape [44j |. 

In conclusion, we note that the theory presented is 
developed for the simplest case of PR generated by the 
uniformly moving point charge. But it may be compara- 
tively easy generalized for radiation of N non-interacting 
point charges, that allows to take into account possible 
coherence effects for all the PR types. Finally, the exter- 
nal field E° may correspond to the field of a multipole 
(multipoles) or even of a non-uniformly moving particle 
(particles), that is of importance, for example, for mod- 
ern laser-plasma based THz radiation sources [t| . 
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Appendix: Exact solution of Maxwell's equations 

In the paper we considered the radiation processes re- 
lated to the far-field part of polarization current density 



field. But along with transverse waves of polarization 
radiation, there may also exist longitudinal waves corre- 
sponding to the zeros of permittivity of the target mate- 
rial. The possibility of existence for such waves is clear at 
least from the fact, that when deriving far- field equality 
([7]) from the general expression (|S|), we used the condi- 
tion r A/ yje(uj). This condition cannot be satisfied 
if e(lu) — > 0, and the ordinary wave zone formulas used 
above become inapplicable. On the other hand, at the 
small enough distances (actually, if r ^ 7A/(27r)) longi- 
tudinal fields may be of importance, even if the value of 
permittivity is not close to zero. Such fields are often 
called wake fields, but the theoretical methods used in 
the corresponding theories noticeably differ from those 
used in the theories of polarization radiation. In particu- 
lar, the physical analogy between the wake fields theory 
and the one of TR and DR was demonstrated relatively 
recently and only for some simplest cases of the ideally 
conducting targets and ultrarelativistic energies of the 
particles (see, for example, @, 

To demonstrate the origin of longitudinal and trans- 
verse waves, it is necessary to derive the general ex- 
pressions for electric field strength and magnetic field 
strength of polarization current density. As we already 
know the exact expression for magnetic field inside the 
medium ([5]), we can find the analogous formula for elec- 
tric field using Maxwell's equations. The simplest way is 
to use the one: 



curlH po '(r,w) + (r , uj) = — j poi (r, u)(A.l) 

c c 



And thus 



W°\v,uj) = ^£-fcurlHP°' - — F° m ), (A.2) 



where = er(w)E°(r, lu) is denoted. As the point of 

observation r is disposed inside the homogeneous medium 
(target of a given shape), e and a are independent of 
coordinates. Then the general expression for magnetic 
field suitable at the arbitrary distances is found from ([5]) 
as 



x (iy/e(u})- 
V c 



r — r 



poi(O) (I 



V T 



dV(A.3) 



Taking curl of this expression according to (|A.2|1 . it is 
possible to derive the following formula for electric field 
strength: 
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uje(uj) 

uj\ 2 3i \Je(w)u) I c 



7 i 5 )} 



i 



r'|2 



r — r / r — r 



r — r \ r — r 



r — r 



r — r 



r — r\ |r — r i 

iy/e(u>)u>\r-r'\/c i 



poi(O) 



r - r' 



(A.4) 



where we used relation ([2]) between conductivity and per- 
mittivity, and integration is performed over the target 
volume. It is the second term in Eq. (|A.4l) 

E- = -i^ T r'(°) = (-^-l)E°, (A.5) 

that describes the longitudinal oscillations. Necessity of 
presence of this term may also be shown via another 
Maxwell's equation: 

divE po ' =47rp po '. (A.6) 

The polarization charge density p po1 is related to the cur- 
rent density with continuity equation: 

pP oi = _L div j P oi = ( divE o + divE poZ ).(A.7) 

ILJ iuj 

Comparing two last equations, we come to the following 
conclusion 

divE"" 1 = Airp po1 = div((^y - l)E ) . (A.8) 

As the integral term in (|A.4|) represents the curl of H po1 , 
divergence of (|A.4I) completely coincides with (|A.8[) . The 
charge density p w = cf° l is sometimes called the den- 
sity of a wake charge [22|]. The field (IA.5|) E™ 1 describes 
the longitudinal waves produced inside the medium by 
the field of the charged particle. Condition e(uj) — > is 
known to be the one of appearance of such oscillations 

num. 

If the condition |r — r' | ^> A/ y/e(ui) is fulfilled, formula 
(IA.4I) transforms into 

c z J \r — r'| L |r — r'| J 

V T 

iWs(u)u\r-r'\/c . -i . 

x . t f3 r / + (- rT -l)E (r jW ). (A.9) 

|r — r| ' 

The field of the uniformly moving charge in vacuum 
E°(r,w) decays exponentially (see e.g. [31(). In an un- 
bounded homogeneous medium, the effective region giv- 
ing a contribution to the integral in (IA.4[) is r e // ~ 
r yX/(2ir). At the distances r ^> -f\/(2ir) the wake field 
may be ignored, and the formula (|A.9|) describes just 
transverse waves of polarization radiation: 

EP^E^-^-exex 

c z r 



x / -f° m exp{i VeR^lr - r'|/c}dV, 

v T 

H poi w H i? = e x E « (A. 10) 

where e = r/r is denoted. As was shown in the 
previous Sections, these simple expressions describe all 
the types of polarization radiation inside the medium. 
The expansion of exponent into r'/r -C 1 series leads 
to the formulas for polarization radiation in the wave 
zone (r ^> 7 2 A/(27r)) or in the so-called pre-wave zone 
(7 2 A/(2?r) > r > 7A/(2?r)). Note that in this paper 
we do not consider the possible interference between the 
fields of polarization currents (|A.3j) . (|A.4[) and the ones 
of the moving particles (formation zone effects). 

In another limiting case of small distances |r — r'| <C 
X/ y / e(uj), or if the value of permittivity is in the vicinity 
of zero: e(w) — > 0, the general formula (|A.4[) describes 
only the longitudinal fields of the oscillating dipoles and 
the wake field: 

C II 

■ppoi ~ -pstat 1 / / r ~ r / r - r .poi(o)\_ 
ws(oj) J V |r-r'| V|r-r'|' J J 

V T 

-r m )y^d^ + (-*- - l)BP(r, W ), (A.11) 

and H po1 ~ jj_ stat = q. The first term here represents 
the static field of the dipole moments with volume density 
p(o) =j(°)/(-ia;) = (e(w)-l)/(47r)E° induced inside the 
medium by the particle's field. Of course, the detailed 
investigation of properties of the longitudinal waves in 
the medium of a finite conductivity (as well as in plasma) 
must include the consideration of the spatial dispersion 
of permittivity. 

It is this field E sta * that acts on the particles being 
in the wake of a leading particle of a bunch moving ac- 
cording to the law r = v(t)t through an obstacle with 
given dielectric properties, or nearby a diffraction radia- 
tion target. If the impact parameter in a problem of DR 
is not too large compared to the wavelength, it is pos- 
sible to consider the particle's energy losses due to the 
force F = eE stat using directly Ea. (IA.ll[) and without 
the application of reciprocity theorem. For example, if 
a bunch moves through the cylindrical channel with ra- 
dius a ^ A in the target of arbitrary conductivity and 
width (see Fig (J), the kick due to the wake fields (|A.11|) 
may be calculated finding the impedances via the stan- 
dard methods. Comparison of the corresponding solution 
for the special case of the ideally-conducting target with 
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well-known result for a waveguide [45J is of interest, but anyway goes beyond the scope of this paper. 
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